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ABSTRACT 


The  present  technical  report  summarizes  the  progress  made  in  considering 
thermal  radiation  in  the  analysis  of  multiphase  turbulent  flow  structure  of  a 
halogren  gas  jet  reacts  with  a  liquid  metal.  The  content  herein  is  based  on 
the  M.S.  dissertation  submitted  by  C.F.  Chern,  Department  of  Mechanical 
Engineering,  University  of  Wisconsin-Milwaukee. 

The  reaction  of  a  halogen  gas  jet  submerged  in  an  akali  liquid  metal  fuel 
forms  a  complex  diffusion  flame.  The  flame  is  turbulent  and  multiphase  in 
nature  and  its  temperature  is  relatively  high.  The  purpose  of  this  study  is 
to  account  for  thermal  radiation  and  the  radiation/turbulent  interaction  in 
analyzing  such  a  diffusion  flame  structure. 

A  method  is  hereby  proposed  to  consider  the  turbulence-radiation 
interaction  in  the  analysis  of  a  multiphase  diffusion  flame.  Application  is 
made  to  study  the  reaction  of  a  halogen  gas  jet  submerged  in  an  alkali  liquid 
fuel  bath.  The  analysis  makes  use  of  a  k-s-g  turbulence  model,  a  local 
homogeneous  two-phase  model,  a  chemical  equilibrium  combustion  model,  a 
probabilistic  approach  and  a  radiation  model  compatible  with  mixture  fraction 
fluctuation.  The  equilibrium  state  relationships  for  instantaneous 
properties,  required  for  mean  scalar  property  predictions  and  structure 
analysis  of  the  reacting  jet,  have  been  generated.  The  energy  equation 
containing  the  radiative  source  term  has  been  solved  simultaneously  with  the 
Favre-averaged  transport  equations.^  Numerical  results  are  presented  for  the 
specific  reactants — chlorine  gas  as  oxidant  and  liquid  sodium  as  fuel.  In 
comparison  to  the  adiabatic  case  in  which  thermal  radiation  is  neglected,  the 
results  reveal  that  thermal  radiation  greatly  affects  the  plume  envelope  and 
the  turbulent  flame  structure. 
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1 .  INTRODUCTION 


The  recent  interest  in  stored  chemical  energy  propulsion  systems  uses  a 
closed  combustor  in  which  thermal  energy  is  generated  by  the  reaction  of  a 
halogen  gas  with  an  alkali  liquid  metal.  The  halogen  gas  is  injected  into  the 
combustor  filled  with  the  alkali  molten  metal,  forming  a  complex,  turbulent, 
multiphase  diffusion-like  flame.  Prior  study  of  this  type  of  flame  is  rather 
limited  (Avery  and  Faeth  1974,  Chen  and  Faeth  1983,  Chan  et  al,  1988)  and  none 
has  yet  considered  thermal  radiation  effects,  in  spite  of  the  high  temperature 
mixture  in  the  flame. 

Prior  studies  on  turbulence-radiation  interaction  in  combustion  are 
confined  to  a  single-phase,  gaseous  flame.  Notable  is  Tamanini's  work  (1977) 
which  assumes  that  the  net  rate  of  volumetric  radiative  heat  loss  is  a  fixed 
fraction  of  the  local  rate  of  the  heat  of  reaction.  The  thermal  radiation 
term  is  in  essence  eliminated  from  the  energy  equation  by  combining  it  with 
the  reaction  heat  source  term  to  circumvent  problems  arising  from 
turbulence-radiation  interactions.  Kounalakis  et  al.  (1988),  recently 
examined  the  effect  of  concentration  fluctuations  on  the  radiation  intensity 
emitted  from  a  gaseous  flame  by  a  stochastic  technique  which  simulates  scalar 
property  variations  along  the  radiation  path.  The  reverse  effect  of  radiation 
on  turbulence  was  not  investigated.  The  present  study  intends  to  account  for 
the  coupled  effects  of  turbulence-radiation  interaction  in  the  analysis  of 
multiphase,  liquid  metal-fuel  combustion. 

Many  existing  radiative  transfer  techniques  (Chan  1987)  cannot  be 
directly  adopted  for  the  analysis  of  a  liquid  metal  combustion  jet  due  to  a 
number  of  difficulties.  8eing  a  turbulent,  multiphase  and  high  temperature 
flame,  the  species  concentrations  in  the  flame,  for  example,  vary  greatly  and 
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fluctuate  rapidly  so  that  only  mean  variables  can  be  derived  from  the  solution 
of  mean  conservation  equations.  Thus  a  compatibility  problem  is  encountered 
when  any  radiative  transfer  technique  is  used  since  radiative  transfer  is 
based  on  instantaneous  variables.  Furthermore,  when  radiative  transfer  is 
considered,  the  flame  becomes  non-adiabatic  and  enthalpy  becomes  a 
nonconserved  scalar.  The  powerful  Shvab-Zel 'dovich  formulation  of  the 
conservation  equations  is  unobtainable  so  that  the  energy  equation  needs  to  be 
solved  explicitly  and  simultaneously  with  other  conservation  equations. 

Another  difficulty  is  associated  with  the  equilibrium  state  relationships 
which  will  become  a  function  of  two  variables  (say  enthalpy  and  mixture 
fraction)  rather  than  one  variable  as  in  the  adiabatic  flame  case.  Thus,  new 
equilibrium  state  relationships  have  to  be  constructed.  Finally,  the 
evaluation  of  mean  scalars  by  a  probability  density  function  (PDF)  technique 
becomes  troublesome  because  the  PDF  is  now  a  function  of  two  variables  as 
well.  All  these  difficulties  will  have  to  be  resolved  in  order  to  include 
thermal  radiation  in  the  analysis  of  a  turbulent  diffusion  flame  in  liquid 
metal  combustion. 

The  primary  purpose  of  the  present  work  is  to  provide  a  methodology  by 
which  the  coupled  effect  of  turbulence-radiation  interaction  in  multiphase 
reacting  flames  can  be  studied.  The  proposed  method  overcomes  the 
difficulties  mentioned  above.  To  illustrate  the  essence  of  the  method  with 
quantitative  results,  a  specific  pair  of  reactants  is  selected — chlorine  gas 
as  an  oxidant  and  liquid  sodium  as  a  fuel.  Since  these  reactants  were 
investigated  previously  (Chen  and  Faeth  1983),  their  results  can  be  used  for 
comparison  in  the  limiting  case  of  an  adiabatic  flame  when  radiative  transfer 
is  neglected. 
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2.  METHOD  OF  ANALYSIS 


2.1  Transport  Equations 

Consider  a  steady,  axisymmetric  submerged  halogen  (or  halogenated)  gas 
jet  discharging  at  a  choked  speed  into  a  large  liquid  alkali  metal  pool.  The 
reaction  between  the  halogen  gas,  the  oxidant,  and  the  liquid  alkali  metal, 
the  fuel,  forms  a  multiphase  turbulent  diffusion-like  flame.  Following  the 
previous  works  (Chen  and  Faeth  1983,  Chan  et  al.  1988),  the  local  homogeneous 
flow  two-phase  model  and  k-c-g  turbulent  model  are  adopted  here.  The 
exchange  coefficients  of  species  are  assumed  to  be  equal  and  thermodynamic 
equilibrium  be  maintained  at  each  point  in  the  flow.  The  Favre-averaged, 
(indicated  by  the  superscript  (~)),  transport  equations  to  be  solved  are 
summarized  in  Table  1.  They  include  the  equations  for  the  mean  axial  velocity 
o',  mixture  fraction  T,  turbulence  kinetic  energy  k,  turbulence  dissipation 
rate  c,  the  square  of  mixture  fraction  fluctuations  g,  and  the  stagnation 
enthalpy  Tf.  Readers  are  refered  to  the  previous  study  (Chan  et  al.  1988)  for 
explanation  of  notations.  The  major  difference  from  the  earlier  studies  lies 
in  the  addition  of  the  energy  equation,  in  which  the  stagnation  enthalpy  is 
the  dependent  variable  and  a  radiation  source  term  has  been  included.  The 
stagnation  enthalpy  can  be  approximated  by  the  static  enthalpy  because  the 
kinetic  energy  is  small  in  comparison  to  the  energy  release  of  the  reaction. 

V«q_  in  the  energy  equation  is  a  Favre-averaged  quantity.  It 
K 

involves  an  integration  of  directional  intensities  over  4ir  solid  angles. 

Each  directional  intensity  can  be  evaluated  stochastically  as  done  by 
Kounalakis  et  al.  (1988).  However,  three-dimensionality  of  the  intensity  and 
the  presence  of  multi-species  in  multiphase  make  this  rigorous  approach 
unpractical  because  of  excessive  computational  time  required  at  each  grid 


-  3  - 


point.  A  great  simplification  can  be  made  if  a  radial  one-dimensionality 
assumption  is  invoked  (Grosshandler  and  Bhattacharjee,  1987),  which  is 
compatible  with  the  current  parabolic,  boundary  layer  type  of  flow.  The 
assumption  is  reasonable  since  flame  is  often  thin  and  elongated  in  shape. 
Further  simplification  can  be  made  when  deemed  necessary  to  reduce  the 
computational  time. 

In  the  present  turbulence/radiation  interaction  (called  coupled) 
analysis,  the  flow  properties  can  be  found  by  solving  the  transport  equations 
for  u,  f,  g,  k,  c  and  "h.  Other  local  mean  scalars  5^,  (such  as  the  mean 
species  mass  fraction,  mixture  temperature  and  density)  can  be  evaluated  by  a 
PDF  averaging  technique  with  some  stipulation.  Since  4>  is  now  a  function  of 
f  and  H  it  is  assumed  (Abou-Ellail  et  al.  1978)  that,  in  as  far  as  the 
evaluation  of  these  mean  scalars  is  concerned,  the  instantaneous  variation  of 
H  and  f  at  the  i-th  grid  point  can  be  approximated  by 


H  =  H  +  f(H.  -  H  )/f .  for  0  <  f  <  f . 

ea  1  49  T  7 


(1) 


=  H.  +  (f  -  f.)(H^  -  H.)/(l-f.)  for  f .  <  f  <  1 

which  is  shown  in  Fig.  1.  Then  the  mean  scalars,  like  Favere  mean  quantities, 
are  evaluated  from 


<0  =  ;  <Kf,H(f,  f^,  H.))  Pp(f,  f.,  g)df  (2) 

0 

where  Pp  is  the  clipped  Gaussian  PDF  profile  (Lockwood  and  Naguib  1975)  and 
()>  is  the  scalar  equilibrium  state  relations  needed  to  be  constructed  next. 
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2.2  Equilibrium  State  Relationships 

Perhaps  the  foremost  complexity  introduced  by  the  inclusion  of  radiative 
transfer  in  the  analysis  is  that  the  mixture  enthalpy  is  no  longer  a  conserved 
property  and  that  other  scalars  at  any  location  will  be  governed  by  two 
variables,  mixture  fraction  f  and  enthalpy  H,  instead  of  f  alone. 

The  equilibrium  mixture  after  reaction  may  be  of  multiphase,  containing 
ideal  or  nonideal  solutions  with  electrolytic  or  non-electrolytic 
constituents.  Its  equilibrium  temperature,  density,  void  fraction  and  mass 
species  fraction  in  all  phases  should  be  calculated  as  a  function  of  the 
mixture  fraction  f  and  the  mixture  enthalpy  H,  over  an  entire  range  of  f  (0  < 
f  <  1)  and  H  (H  <  H  <  H  ).  The  subscripts  o  and  »  designate  the 
nozzle  exit  (f=l)  and  bath  fluid  (f  =  o)  conditions.  These  equilibrium  state 
relationships  can  be  constructed  using  the  CEC-NMS  equilibrium  code  recently 
developed  in  this  laboratory  (Chan  et  al.  1988). 

2.3  Summary  of  the  Methodology 

Starting  from  the  nozzle  inlet  condition,  the  Favre  mean  radiative  source 

term  a^  the  first  axial,  downstream  location  is  calculated  by  a  stochastic 

technique  based  on  upstream  radial  conditions.  This  is  calculable  in  view  of 

the  radial  one-dimensionality  assumption.  With  the  resulting  radiative  source 

term,  the  transport  equations  listed  in  Table  1  are  solved  to  yield  dependent 

variables,  like  T.  and  H!. ,  etc.,  at  the  first  downstream  location.  With 
1  i 

equations  (1-2)  and  the  equilibrium  state  relationships  <j>  =  4>  (f,H),  other 
mean  scalars  not  yet  solved  are  calculated.  If  desired,  iteration  can  be  made 
using  new  p  and  7«q  based  on  the  downstream  condition.  The  procedure  is 
repeated  at  the  next  downstream  location.  Since  the  radiative  source  is 
calculated  by  a  stochastic  technique,  the  effect  of  turbulence  on  radiation  is 
accounted  for.  More  importantly,  when  the  radiative  source  term  is 
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subitituted  into  the  energy  equation  which  is  then  solved  simultaneously  with 
the  turbulence  momentum,  K  and  c  equations,  they  are  coupled  through  u,  v, 
and  p  present  in  those  equations.  The  effect  of  radiation  (through 
p  in  u,  k  and  c  eqs.)  on  turbulence  and  the  effect  of  turbulence  (through 
u,  V  and  in  f  and  H  eqs.)  on  radiation  are  further  considered.  The 
latter  is  apparent  since  local  and  H.  are  used  (see  Eq.  1)  to  evaluate 
the  fluctuating  f  and  H  quantities,  which  in  turns  are  employed  to  calculate 
the  instantaneous  species  compositions,  temperature  and,  consequently,  the 
instantaneous  radiation  properties  and  intensity.  Finally,  phase  change  and 
chemical  reaction  are  coupled  to  radiation  and  turbulence  through  the  use  of 
the  equilibrium  state  relationships.  Thus  the  method  proposed  indeed  accounts 
for  turbulence,  radiation,  phase  change  and  reaction  interactions. 

3.  CHLORINE  JET  SUBMERGED  IN  LIQUID  SODIUM 

To  illustrate  the  essence  of  the  methodology,  consider  the  steady 
reaction  of  an  axisymmetric  chlorine  gas  jet  discharged  at  a  choked  speed  from 
below  into  a  large  molten  sodium  pool.  The  needed  information  for 
computation  is  presented  next. 

3.1  Equilibrium  State  Relationships 

For  Cl^/Na  reactions,  the  possible  species  in  the  gas  phase  are 

‘^^2(G)'  ’^^G)'  '^®2(G)'  ^®2^^2(G)-  Similarly, 

the  species  in  the  liquid  phase  are  and  NaCl^j^^,  and  the  species  in 

the  solid  phase  are  Na^^j  and  NaCl^^j.  The  new  equilibrium  state 
relationships  of  the  mixture  temperature,  density  and  mass  fractions  of  all 
species  have  been  constructed  as  a  function  of  f  and  H,  and  are  illustrated  in 
Figs  2  to  19.  The  conditions  before  the  reaction  are  T^  =  298  K  (the 
oxidant  temperature),  T  =  1130  K  (the  fuel  temperature)  and  =  1 

QD  0 
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in  the  flame.  Then  on  the  basis  of  mixture  fraction  profiles,  the  flame 
envelope,  as  shown  in  Fig.  23  and  Table  2,  is  subdivided  into  two  segments  in 
the  axial  direction  and  five  segments  in  the  radial  direction  for  a  total  of 
ten  sublayers.  The  mean  radiation  properties  in  each  sublayer  are  estimated 
next. 

In  the  Na/Cl^  reaction,  the  primary  radiation  participating  species  are 
NaCl^j^j,  NaCl^jj,  NaCl^gj  and  Na^Cl^^gj.  NaCl^^j  mass 
fraction  is  relatively  small  and  can  be  neglected. 

To  calculate  the  absorption  and  scattering  coefficients  attributable  to 

the  liquid  phase,  the  mass  of  each  liquid  species  and  NaCl^l^^)  in 

each  sublayer  is  calculated  as  presented  in  Table  3.  The  liquid  species  are 

assumed  to  exist  in  droplet  form.  Since  the  droplet  diameter  and  density 

can't  be  predicted  accurately,  the  diameter  is  estimated  by  a  Weber  number 

2 

droplet  break-up  criterion  With  the  droplet 

size  and  liquid  mass  determined,  the  droplet  density  of  species  i,  N.,  in 
each  sublayer  can  be  calculated  (see  Table  4).  The  absorption  and  scattering 
coefficients  are  then  evaluated  using  the  large  particle  limit  model  (i.e. 
irO./x  >  5,  X  being  radiation  wave-length). 


'hi 


ir  o;/4,  K  =  N. 
1  s.  1 


^hi 


v  D^/4 


(3) 


For  a  liquid  sodium  at  1130  K,  c.  ,,  .  and  p.  ,,  ...  are  found  from 

h,Na(L)  h,Na(L) 

Chan  and  Tseng  (1979).  As  to  a  NaCl  liquid  droplets,  its  value  is 
approximated  by  that  of  a  large  ideal  dielectric  sphere.  The  absorption  and 
scattering  coefficients  of  the  liquid  droplets  are  summarized  in  Table  5. 

The  absorption  coefficient  of  a  gas  species  is  calculated  from  the 
sublayer  mean  beam  length  and  the  species  emissivity  by 
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(4) 


In  the  absence  of  data  on  emissivities  of  NaCl,^.  and  Na»Cl-.-.  their 

(6)  2  2(G) 

orders  of  magnitude  are  arbitrarily  approximated  by  that  of  a  triatomic  gas, 

say,  water  vapor,  at  the  same  temperature  and  pressure  path  length  (P  .L  ) 

ei  e 

conditions.  The  estimated  absorption  coefficient  of  the  gas  phase  is  given  in 
Table  6. 


Finally,  by  summing  up  the  contributions  from  all  species  in  all  phases, 
the  mixture  absorption  and  scattering  coefficients  in  each  of  the  ten 
sublayers  are  yielded  and  summarized  in  Table  7.  The  scattering  coefficient 
is  included  here  for  later  study. 

Similarly,  a  simplified  radiative  source  expression  is  desirable.  Under 
the  planar  approximation,  the  source  term  in  an  absorbing,  emitting 

and  scattering  medium  is  [Sparrow  and  Cess,  1978]: 


4K  K 

^•‘»rad  =  «b(^)  -  -  2^0,1  '  2^0,2^  ^^a  " 


'o,l 


^0,2 


G(X)  = 


ciepi  (1  -  c^)c2%2 

1  -  (1  -  c^)(l  -  Cg) 


ce  +(l~c)ce 
2  b2  ^  2^  1  bl 

1  -  (1  -  c^)(l  -  c^) 


(5) 


(6) 


(7) 


2ir  J  i'^(0,p)e“’'''’'dw  +  2ir  J  i'e^’^o 
0  0 


°  '^a  S 

_« _ /  4.  \  .  _  >_ 


(8) 


a  s' 


For  the  optically  thin  limit,  the  final  form  of  Eq.  (5)  becomes 


=  2K  [2e  -  (q„  ,  +  q„  J] 
rad  a  b  o,l  o,2 


(9) 


If  the  wall  emissivities  (cl  and  c2)  and  temperature  and 

the  same,  then  -  ^bl  '  ®b2  '  'b.surr  ='’‘’ 

’•Vad  -  «a'’<T''  -  0 


(10) 


in  which  surface  properties  and  the  scattering  coefficient  vanish  nicely.  In 
Eq.  (10),  is  the  bath  temperature  and  K^,  the  absorption 

coefficient,  has  been  calculated  and  listed  in  Table  7.  Equation  (10),  though 
derived  from  the  planar  layer,  is  equally  applicable  for  the  cylindrical 
geometry.  This  approximation  is  adopted  here  since  the  current  axisymmetric 
flame  is  elongated  and  optically  thin  radially.  Gross  estimation  of  K 

d 

values  seems  to  render  unnecessary  the  evaluation  of  the  Favre  averaging  of 

(V»q  )  by  a  rigorous  stochastic  technique.  Rather  the  local  K  and  T 
K  a 

are  approximated  by  the  mean  values  of  the  local  sublayer  to  reduce  the 
excessive  computation  time  in  this  preliminary  investigation. 


3.3  Boundary  and  Initial  Conditions 

The  uniform  Cl^  inlet  conditions  at  the  nozzle  (diameter,  d  =  7.9  x 
10'%  exit  are  =  2.898  kg/m^,  T^  =  298  K,  “u^  =  203.4m/s,  =  (O.OBu^)^, 

=  (0.03tr^)^/d,?g  =  1,  =  0,  'Hq(T^)  =  64.22  J/kg.  The  sodium  bath 

5 

conditions  are  T^  =  1130  K  and  =  10  Pa.  In  the  outer  free  boundary  (the 

Nth  grid  point)  of  the  reacting  plume,  'uj^  =  0.0005  “Ug,  fj^  =0.0005,  gj^  =  0, 

^  6 

H|^(T^)  =  1.188  X  10  J/kg  while  kj^  and  are  obtained  by  solving  the 
governing  equations  while  neglecting  the  radial  convective  and  diffusive  terms 
as  follows: 
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^d.N  =  ^.N  -  “^u.N 


‘d,N  “  %,N  ^g2‘'u.N  “u.N^ 


(11) 

(12) 


/'✓ 


In  the  above,  H  and  H.,  are  enthalpies  of  pure  chlorine  and  pure  liquid 
sodium,  subscripts  (d,N)  and  (u,N)  denote  downstream  and  upstream  of  Nth  grid 
point,  and  ax  is  the  step  size  in  the  jet  axial  direction. 

The  gradient  of  the  dependent  variable  is  set  to  zero  for  the  inner 
symmetrical  boundary  along  the  center  jet  axis,  d%/dr  =  0,  where  ♦  =  V, 
f,  g,  k,  c  and  H  from  the  adiabatic  flame  structure  solution. 


4.  METHOD  OF  SOLUTION  AND  RESULTS 


4.1  Adiabatic  Flame 

The  flame  structure  computations  were  performed  using  a  version  of  GENMIX 
code  (Spalding,  1977)  modified  by  our  team.  Before  computing  the  coupled 
radiation/turbulence  interaction,  the  adiabatic  flame  structure  was  computed 
first,  partly  for  radiation  property  evaluation  and  partly  for  comparison  with 
the  existing  solution  (Chen  and  Faeth  1983).  The  variations  of  the  mean 
properties  of  temperature,  nondimensional  velocity,  mixture  fraction,  void 
fraction  and  the  overall  mass  fraction  of  each  phase  along  the  dimnesionless 
axis  location,  x/d,  are  shown  in  Fig.  24.  The  highest  temperature  along  the 
axis  is  about  2950K  at  x/d=26.  This  is  lower  than  the  highest  possible 
instantaneous  temperature  because  of  the  use  of  a  clipped  Gaussian  probability 
density  function.  According  to  the  axial  void  fraction  profile  in  Fig.  24, 
the  flame  length  (the  total  penetration  length)  extends  to  x/d=150.  By 
referring  to  Figs.  24  and  25,  it  is  seen  that  reaction  is  completed  at  the 
position  where  the  temperature  reaches  its  maximum  which  is  nearly  at  the 
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stoichiometric  condition.  The  axial  variations  of  the  mean  species  mass 
fraction  in  different  phases  are  also  plotted  in  Figs.  25  to  27.  In  the  gas 
phase,  most  of  species  are  Cl^.  Na  and  NaCl.  NaCl  has  the  largest  amount  at 
x/d=26  where  temperature  reaches  the  maximum.  This  is  understandable  since 
the  reaction  is  almost  a  simple  and  fast  type  of  reaction  at  the 
stoichiometric  condition.  Most  of  the  liquid  phase  appear  when  x/d  is  greater 
than  25  in  which  temperature  is  lower  than  the  boiling  temperatures.  As 
mentioned  before,  the  solid  phase  occurs  only  when  the  mixture  fraction 
approaches  1  where  the  mixture  temperature  is  much  less  than  the  melting 
temperature  of  NaCl  (1073K).  This  is  shown  in  Fig.  27.  Physically,  solid 
phase  exists  only  near  the  nozzle  location  where  mixture  fraction  approaches  1 
and  temperature  is  low.  The  potential  core  where  velocity  is  approximately 
equal  to  or  greater  than  99  percent  of  initial  velocity  terminates  at  x/d=12. 
All  the  calculations  are  in  good  agreement  with  that  of  Chen  and  Faeth  [1983]. 

The  radial  variations  of  the  mean  properties  of  temperature, 
nondimens ional  velocity  and  mixture  fraction  are  also  presented  in  Figs. 

28-31.  The  profiles  of  the  mean  mixture  fraction  and  velocity  are  nearly 
identical  because  buoyance  effect  is  neglected  in  the  source  term  of  momentum 
equation.  The  form  of  the  momentum  equation  without  a  source  term  is  similar 
to  that  of  mixture  fraction  equation.  The  slight  difference  between  mean 
mixture  fraction  and  velocity  is  due  to  different  effective  viscosity 
and  boundary  conditions  used  in  both  equations. 

Additional  calculations  for  a  pure  chlorine  injected  into  a  pure  liquid 
sodium  at  atmospheric  pressure  indicate  that  the  flow  structure  and  scalar 
properties  are  significantly  controlled  by  the  ambient  temperature  because  of 
subcooling  effect.  Decreasing  liquid  temperature  will  increase  condensation 
which  in  turn  reduces  the  flame  length  as  defined  by  the  void  fraction 
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profile.  The  effect  of  subcooling  is  shown  in  Figs.  32-33  and  Table  8.  In 

Fig.  22,  T  is  chosen  as  902K  while  in  Fig.  33,  T  is  chosen  as  674K. 

00  00 

T„,  P  and  other  conditions  are  the  same  as  the  case  of  T  =1130K. 

0  o»  00 

From  the  above  Figures  and  Table,  we  can  see  the  significant  effect  of 
subcooling  on  turbulent  diffusion  flames  of  a  gaseous  oxidizer  injected  into  a 
liquid  metal  fuel.  The  total  penetration  length  will  extend  to  infinity  as 
the  liquid  sodium  approaches  its  boiling  temperature  because  of  the 
diminishing  driving  potential  for  condensation  (i.e.,  no  subcooling  effect). 

4.2  Nonadiabatic  Flame  with  Radiation 

The  computational  results  of  the  coupled  turbulence/radiation  interaction 
analysis  are  given  in  Fig.  34-37  which  show  that  the  flow  structure  deviates 
from  that  of  adiabatic  case  presented  above.  In  Fig.  34,  the  total 
penetration  length  of  the  plume  (defined  by  the  length  when  the  mixture 
fraction  reduces  to  0.005  which  approximately  corresponds  to  zero  void 
fraction)  extends  to  x/d=123,  which  is  about  18  percent  shorter  than  the 
adiabatic  one.  The  highest  temperature  along  the  center  axis  is  about  2810  K 
at  x/d=23.  This  is  lower  than  the  adiabatic  case  of  2950  K.  The  profiles  of 
mixture  fraction  f  and  velocity  u/u^  remain  similar  even  when  radiation  is 
accounted  for  in  turbulent  reacting  flow.  There  exists  only  the  gaseous  phase 
between  x/d=16  and  x/d=25  where  the  mixture  temperature  exceeds  the  boiling 
temperature  of  Na  and  NaCl.  In  Figs.  34-37,  all  quantities  shift  to  the  left 
due  to  the  effect  of  radiation/turbulence  interaction.  The  differences  in  the 
radial  temperature  and  velocity  profiles  at  different  x/d  between  coupled  and 
adiabatic  analyses  are  presented  in  Figs.  38-39.  Figure  38  shows  that  the 
velocity  difference  is  small  at  small  x/d  locations.  This  is  because  the 
radiation  effect  is  insignificant  near  the  injector  where  the  jet  stream 
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temperature  is  still  low.  The  difference  becomes  progressively  larger  away 
from  the  injector.  Close  to  the  injector,  due  to  the  radiative  heating  by  the 
surrounding  warmer  liquid  sodium,  the  jet  stream  along  the  jet  axis  is  heated 
up  above  the  adiabatic  case  as  shown  in  Fig.  39(a)  and  (b).  But  then  as  the 
jet  temperature  increases,  the  radiative  transfer  from  the  jet  to  the 
surrounding  becomes  important  and  the  mixture  jet  temperature  eventually  drops 
below  the  adiabatic  one  as  illustrated  in  Figs.  39(c)  and  (d).  By  referring 
to  Figs.  24,  34  and  39(c),  the  maximum  temperature  difference  between  the 
coupled  case  and  the  adiabatic  case  is  about  300  K  at  x/d=25.  In  Figs.  39(c) 
and  (d),  due  to  radiation,  the  temperature  of  the  coupled  case  is  lower  than 
the  adiabatic  case  near  the  center  axis.  Also  shown  in  Figs.  39(c)  and  (d), 
the  location  of  second  grid  point  y^  of  te  coupled  case  is  slightly  farther 
away  from  the  center  axis  than  that  of  the  adiabatic  case.  This  is  because  in 
th  6ENMIX  formulattion  [Spalding,  1977,  pp  112-113],  the  y^  location  can  be 
readily  shown  as 


.!VV!i1.i/2 

2  ■  ^  (pu)^^  > 


(11) 


where  =  1  in  the  present  axisymmetric  flow,  is  a  pre-determined 
dimensionless  stream  function  at  location  2  and  (¥^-fj)  is  the 
entrainment  rate.  For  the  coupled  case,  since  the  entrainment  rate  is  larger 
and  the  velocity  u^  is  smaller  than  the  adiabatic  case,  y^  is  larger 
accordingly.  In  addition,  all  quantities  at  y^  in  the  GENMIX  formulation 
are  inherently  identical  to  those  at  y^  (=0).  This  explains  the  flatness  of 
the  curves  shown  in  Figs.  38  and  39. 

The  differences  of  the  jet  profile  and  plume  profile  are  shown  in  Figs. 
40  and  41,  respectively.  The  jet  profile  is  defined  as  the  envelope 
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connecting  all  the  Nth  grid  points  (the  furthest  grid  point  from  the  jet  axis) 
while  the  plume  profile  is  defined  similar  to  the  total  plume  length  as  the 
locations  where  f  is  reduced  to  0.005.  At  such  low  values  of  f,  the  void 
fraction  is  practically  zero.  Figure  40  shows  that  the  jet  profile  is  not 
greatly  affected  by  the  radiation/turbulence  interaction.  In  Figure  41,  the 
plume  profile  of  the  coupled  case  begins  to  shrink  around  x/d=80,  much  shorter 
than  the  adiabatic  case  which  occurs  around  x/d=130.  Due  to  the  collapsing 
plume  and  therefore,  the  void  fraction,  the  overall  jet  profiles  also  begin  to 
shrink  to  some  extent  at  those  x/d  values  as  shown  in  Fig.  40.  This 
phenomenon  is  attributed  to  radiative  heat  loss  which  causes  a  lower  mixture 
temperature.  Decreasing  mixture  temperature  will  reduce  mixture  volume  under 
constant  pressure  condition  and  therefore  shrink  the  jet  envelope. 

5.  CONCLUSION 

A  method  has  been  proposed  to  consider  thermal  radiation  in  the  analysis 
of  a  turbulent  multiphase  diffusion  flame.  A  specific  example  is  presented 
for  a  reacting  chlorine  gas  jet  submerged  in  a  molten  sodium  pool.  In  the  ab¬ 
sence  of  data,  thermal  radiation  properties  of  Na/Cl^  reaction  products  were 
estimated.  New  equilibrium  state  relationships  for  scalar  properties  were 
presented  as  a  function  of  enthalpy  and  mixture  fraction.  Results  of  flame 
structure  including  flame  envelope  have  also  been  presented.  Significant 
differences  in  flame  structures  were  found  between  the  adiabatic  flame  and  the 
non-adiabatic  flame  which  accounts  for  radiative  heat  losses.  The  results  of 
the  computation  shows  a  strong  turbulence/radiation  interaction  even  under 
the  optically  thin  condition. 
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TABLE  1.  Governing  Equations 


VI  cl  g2  gl  gZ  k  c  r  g  H  c 

0.09  1.44  1.84  2.8  1.84  1.0  1.3  0.7  0.7  0.7  0.7 

.  =  pC  k^/c 
t  p 


locations  of  regions  in  radial  end  axinl  directions 


The  number  in  the  parenthesis  corresponds  to  each  mixture  fraction 
profile  and  is  numbered  outward  radially 


sublayerNo 

1 

3 

4 

5 

6 

7 

8 

9 

lOj 

Nai 

179 

0 

0 

0 

0 

684 

0 

0 

0 

u 

NaCh 

151 

7.26 

0 

0 

0 

558 

24.2 

0 

0 

El 

NaCla  +  Na^CliG 

11.9 

74.2 

96.2 

60.5 

2.86 

41.6 

238 

264 

110 

El 

unit :  10~^p 


Table  3:  The  mass  of  each  sublayer 


sublayerNo 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Nai^ 

19124 

0 

0 

0 

0 

73077 

0 

0 

0 

0 

NaCh 

4247 

268 

0 

0 

0 

20590 

593 

0 

0 

0 

Table  4:  The  number  of  particle  in  each  sublayer 


sublayerNo 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Ka,Sa,L 

0.00036 

0 

0 

0 

0 

0.00133 

0 

0 

0 

0 

0.027 

0 

0 

0 

0 

0.101 

0 

0 

0 

0 

^aJ^aCl,L 

0.006 

0.0003836 

0 

0 

0 

0.029 

0.000849 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

unit  of  Ka  and  A’,:cm“‘ 

Table  5;  The  absorption  and  scattering  coefficients  of  liquid 


unit :  10”®cm“' 


Table  6  :  The  absorption  coefficient  of  gas 


Sublayer 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

XjD 

0-35 

0-35 

0-35 

0-35 

0-35 

35-70 

35-70 

35-70 

35-70 

35-70 

f 

0-.3 

.3-.45 

.45-.7 

.7-.95 

.95-1 

0-.3 

.3-.45 

.45-.7 

.7-.95 

.95-1. 

Ka 

26.8 

24 

14 

12.6 

59.8 

39.3 

11.5 

6.8 

6.6 

108.4 

K 

27 

0 

0 

0 

0 

101 

0 

0 

0 

0 

Table  7 ;  The  mixture  absortion  and  scattering  coefficients  in  combustion  field 


1130 

‘902 J 

674 

FlameLength,  x/d 

150 

25 

10.2 

Table  8:  Comparison  of  subcooling  effect 
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Figure  lo:  Y(NA,G)  distribution  with  different  f  and  H 


Figure  11:Y(NA2,G)  distribution  with  different  distribution  f  and  H 


Figure  14:  Y(NA,L-1)  distribution  with  different  f  and  H 


Figure  15:  Y(NACL,L-1)  distribution  with  different  f  and  H 
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Figure  18:Y(NA,S)  distribution  with  different  f  and  H 


Figure  19;Y(NACL,S)  distribution  with  different  f  and  H 
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Figure  22;  Mass  fraction  distribution  with  mixture  fraction 
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Figure24;  Variations  of  mean  properties  of  temperature,  mixture  fraction,  velocity, 
void  fraction,  gas  phase,  liquid  phase  and  solid  phase  along  the  central  axis 
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Figure  26:  Variations  of  mean  properties  in  liquid  phase  along  the  central  axis 
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Figure  28; Variations  of  temperature,  velocity  and  mixture  fraction  at  x/d=b 
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Figure  30;  Variations  of  temperature,  velocity  and  mixture  fraction  at  x/d=15 
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Figure  32:  Variations  of  mean  properties  along  the  central  axis  due  to  subcooling 
effect 
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Figure  33;  Variations  of  mean  properties  along  the  central  axis  due  to  subcooling 
effect 
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Figure  34:  Variations  of  mean  properties  of  temperature,  mixture  fraction,  velocity, 
void  fraction,  gas  phase,  liquid  phase  and  solid  phase  along  the  central  axis  with 
radiation 
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Figure  40 :  Comparison  of  jet  profile  between  coupled  and  radiation  analyses 
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Figure  41;  Comparison  of  plume  length  between  adiabatic  and  coupled  analyt 
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